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Abstract. We investigate temperature and entropy profiles of 13 nearby cooling flow clusters observed with the EPIC cameras 
of XMM-Newton. When normalized and scaled by the virial radius the temperature profiles turn out to be remarkably similar. 
At large radii the temperature profiles show a clear decline starting from a break radius at ~ 0.1 rvir- The temperature decreases 
by ~ 30% between 0.1 rvir and 0.5 rvir. As expected for systems where non-gravitational processes are of great importance, 
the scale length characterizing the central temperature drop is not found to be proportional to the virial radius of the system. The 
entropy of the plasma increases monotonically moving outwards almost proportional to the radius and the central entropy level 
is tightly correlated with the core radius of the X-ray emission. The dispersion in the entropy profiles is smaller if the empirical 
relation S (X T«'«5 is used instead of the standard self-similar relation S (x T and, as expected for cooling flow clusters, no 
entropy cores are observed. 

Key words. Galaxies: clusters: general - cooling flows - X-rays: galaxies:clusters 



Oh; 

6 



X 



1. Introduction 

Temperature profiles of the X-ray-emitting intracluster 
medium (ICM) have several important implications. Derived 
from X-ray data together with the gas density distribution, they 
are critical in deriving the total gravitational mass of galaxy 
clusters, which combined with the estimated gas mass and pri- 
mordial nucleosynthesis calculations, provides constraints on 
the cosmological density parameter (e.g., [White etal. (1993)^ . 
Temperature profiles provide useful knowledge to comple- 
mentary observations such as the Sunyaev-Zel'dovich effect 
and affect the determination of cosmological parameters like 
the Hubble constant (e.g.,'Rephaeli (1995)*. They also provide 
valuable information on the thermodynamic history of galaxy 
clusters, since they are fundamental in determining the ICM 
entropy distribution, and the latter has been shown to be a 
very powerful tool to study non-gravitational processes such 
as radiative cooling, preheating and feedback from supernovae 
and active galactic nuclei (e.g., |Voit (2 004) l. 
Various observational studies have found different and con- 
flicting results regarding temperature gradients in the outer 
regions of galaxy clusters. Markevit ch et al. (1 998) analyzed 
projected radial temperature profiles for 30 clusters observed 
with ASCA, finding that nearly all the clusters in their sample, 
including cooling flow clusters, show a significant decline at 



large rad iijlmii^^^l 999)| pointed out that the claim by 
|Markevifche^^0998)r of a declining temperature profile 
is in contrast with the findings of other authors even for 
the same object and strengthened the claim of isothermality 
by investigating ROSAT PSPC color profiles. White (2000l 
studied a sample of 106 clusters observed with ASCA and 
found that 90 percent of the temperature profiles are consistent 
with isothermality. [Irwin & Bregman (2000)[ derived radial 
temperature profiles for 1 1 clusters observed with BeppoSAX 
and found generally flat or slightly increasing temperature 
profiles, therefore ruling out a declining temperature profile. 
On the other hand Finoguenovet al. (2001b)^ found dec lining 
profiles similar to those of [MarkevitdTet al. (1998)[ from 
the analysis of ASCA/SIS data of 18 relaxed cool clusters. 
| De Grandi & Molendi (2002), analyzed the projected tempera- 
ture profiles for 1 1 cooling flow and 10 non-cooling flow clus- 
ters observed with BeppoSAX. They find a rapid decline of the 
profiles beyond ^ 0.2rigo and that non-cooling flow clusters 
appear to be well described by an isothermal core extending 
out to ^ 0.2ri8o, while the well known signature of cooling 
flow clusters is the temperature decline towards the center. 
Despite the agreement between [Markevitch et al. (1998) 



'Finoguenov et al. (2001b")' and De Grandi & Molendi (2002) 
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the behavior of cluster temperature profiles beyond ~ 0.2ri8o 
remains controversial, mainly because even with large samples 
the number of clusters where the temperature is derived 
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beyond O.Srigo is very small. 

In the early phases of the XMM-Newton mission, 
[Arnaud et al. (200131 found that the temperature profile 
of A 1795 is essentially isothermal up to ~ 0.4r2oo- 
[Pratt et al. (200 1)| find a flat temperature profile up to 
~ 0.5r2oo for A 2163 and an indication of a mild decline at 
larger radii. XMM-Newton observation of the distant cluster 
RX J1120.1+4318 indicate a flat profile up to O.Srvir 
l |Arnaud et al. (2002)| l, and no sharp temperature drop is 
found in A 1835 up to ~ 0.7rvir ( |Majerowicz et al. (2002^ . 
[Pratt & Arnaud (2002)] find a declining temperature profile in 
A 1413. The gradient is modest compared to the results of 
IMarkevitch et al. (1998)| and |De Grandi & Molendi (2002)] 
the temperature decreases by ~ 20% between 0.1r2oo and 
0.5r2oo- The temperature profile of A 1983 derived within 
^ 0.4r2oo by 'Pratt & Arnaud ( 2003)|does not s how a decline 
and for A 1650 Takahashi & Ya mashita (2003)| find a temper- 
ature decline of ^ 20% between 0.1 r2oo and 0.4r2oo- Given 
these controversial results on temperature gradients at large 
radii, any further insight is welcome. 

In this paper we thus investigate radial temperature profiles 
of nearby cooling flow clusters using the results from the 
spatially resolved spectra taken with the EPIC cameras of 
XMM-Newton presented in Ka astra et al. (2004)1 We present 
an accurate description of temperature and entropy profiles at 
all the observed radii and investigate their scaling properties. 
We additionally provide information on the "uncorrected" 
luminosity-temperature relation for cooling flow clusters. 
The layout of this paper is as follows. In Sect. |2] we briefly 
summarize the procedure and results which are essential for 
our analysis (see [Kaastra et al. (20(j4)| for the description of 
the sample and the extensive presentation of the data analysis) 
and discuss the selection of temperature bins. In Sect. |3l we 
compute, for each cluster, gas and total gravitational masses, 
and the virial radius needed for the scaling of temperature 
profiles presented in Sect. |3 along with the comparison with 
previous work and numerical simulations. In Sect.|5lwe focus 
on the scaling of entropy profiles and present the luminosity 
temperature relation derived from our sample. In Sect. |6l 
we summarize the results of our analysis and present the 
conclusions. 

Throughout this work we use, depending on the con- 
text, a SCDM (fi,n = 1.0, r^A = 0.0) cosmology with 
Ha = 50 kms-^Mpc-i and/or a ACDM (n^ = 0.3, 
JIa = 0.7) cosmology with Hq = 70 kms^^Mpc^^. These 
will be denoted SCDM50 and ACDM70, respectively. The 
main reason for presenting results for both cosmologies is 
that in the literature often results for the SCDM50 cosmology 
are listed, while only recently more people start using the 
ACDM70 model. This makes a comparison of our results with 
other work more easy. 

2. Sample 

The sample consists of 17 clusters with cooling flows listed 
in Table [2 Kaastra et al. (2004)| derived deprojected (3D) and 
projected (2D) radial temperature profiles for all the clus- 
ters in this sample. From the analysis of deprojected spectra 



[Kaastra et al. (2004)[ give the results of single-temperature and 
multi-temperature fittings. Even though the single-temperature 
model does not provide the best fits to the inner cooling parts in 
several clusters, we will use the results from the single-phase 
model for the following reason. The multi-temperature fitting 
indicates that the gas in the cooling flow region is multiphase 
and at the same time that the emission measure distribution 
is dominated by the emission from the hottest gas at each ra- 
dius (see Kaas tra et al. (2004)| Table 7). Using the results of 
the multi-temperature fitting and assuming pressure equilib- 
rium between the gas phases we compute filling factors and 
total gas density for each shell. We then compute the effective 
temperature of the gas mixture in each shell and find that, for 
all clusters and all shells, it agrees within less then a few per- 
cents with the shell temperature from the single-temperature 
modelling. 

Spectral fitting is done over the full 0.2-10 keV range. The 
energy range below 0.2 keV is currently too poorly calibrated 
to be useful for spectral analysis, while above 10 keV most 
of our spectra lack sufficient flux. In general, temperatures 
and gas densities are computed for the innermost 8 shells (see 
[Kaastra et al. (2004), Table 3 for the boundaries between the 
shells), except when the data become too noisy in the outer- 
most shell. 

In the derivation of the total gravitational mass from the 3D 
temperature and density profiles and for the study of temper- 
ature gradients at large radii, it is crucial to select bins with a 
robust temperature estimate. For Virgo and Perseus the inner- 
most shell is excised because of the presence of a strong X-ray 
emitting AGN. At high energies, the background as measured 
with XMM-Newton is dominated by a power law contribution 
caused by soft protons. This component is time variable and 
we estimate that we can determine it with a maximum uncer- 
tainty of 10% for each observation. As in the outer regions of 
the cluster the X-ray flux drops rapidly, this soft proton com- 
ponent dominates (see 'Kaast ra et al. (20 04) Sect. 2.3 for a de- 
tailed presentation on the background subtraction). For this rea- 
son we exclude temperature bins in the outermost parts of the 
cluster as follows. Considering that the background uncertainty 
is of the order of 10%, a bin with temperature kT is selected 
if 10% of the net flux is larger or equal to 10% of the back- 
ground in the energy band containing the turnover region of 
the bremsstrahlung spectrum, i.e. the one closer to the energy 
interval kT - 2kT. Concerning the choice of the energy band in 
the selection procedure, we notice that particular carefulness is 
needed in the selection of the outermost bins with low temper- 
atures. We have applied our selection criterion using different 
energy bands close to the energy interval kT - 2kT and find that 
if flux at too low energy compared to kT is taken into account, 
bins that would otherwise be excluded, satisfy the selection cri- 
terion. This effect can be attributed to the soft excess, since it 
is affecting the outermost bins only and it is strong in clusters 
known to show soft excess (A 2052, Sersic 159—3, MKW 9 and 
MKW 3s). A detailed modelling of the soft excess is beyond the 
scope of this work and we prefer to exclude bins without robust 
temperature estimates even though their inclusion would allow 
us to trace temperature profiles at much larger radii. The se- 
lection is performed using both pn and MOSl/2 data. We find 
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the selection with pn to be more conservative than the one per- 
formed on MOSl/2 data and therefore use the resuhs of the 
former. In Tabled we Hst i?out> the value of the radius at the 
outer end of the last radial bin considered after the selection. 
The inputs of our study are therefore: gas bolometric luminos- 
ity, gas density and temperature in each shell within i?out- 



3. Gravitational masses 

In order to make a fair comparison between clusters, the var- 
ious quantities such as gas temperature and entropy must be 
evaluated at constant fractions of the virial radius, rather than 
at fixed radii imposed by data limits. We therefore estimate the 
total mass assuming spherical symmetry and hydrostatic equi- 
librium and compute r^, the radius within which the mean in- 
terior density is A times the critical value. 
We estimate the total gravitational mass by making direct use 
of the 3D gas temperature and gas density from the single phase 
analysis. We invert the equation of hydrostatic equilibrium, i.e. 



Mtot(< r) 



kT{r) /d\np{r) 
^inipG V d\nr 



dlnT(r) 
dlnr 



(1) 



where G is the gravitational constant, k the Boltzmann con- 
stant, p the gas density, T its temperature and /imp is the mean 
particle mass of the gas (we assume ^ ~ 0.61). Since the mass 
contribution from galaxies is small we neglect it, and there- 
fore Aftot(< r) = A/dm(< r) + Mgas(< r), i.e. the total 
gravitational mass within a sphere of radius r is given by gas 
plus dark matter mass. Using the three-dimensional gas density, 
we select the dark matter mass model that reproduces better 
the deprojected temperature profile. As a cosmologically moti- 
vated dark matter mass model, we consider the integrated NFW 
( [Navarro et al. (1996)) dark matter profile: 

MuM{<r)=iTTrlp,- \_ , ^ , V'^ , (2) 



ln(l 



c/(l 



where pc is the critical density. The scale radius rg and the con- 
centration parameter c are the free parameters. 
We make use of the NFW dark matter density profile because 
the choice of an analytical function allows us to extrapolate the 
dark matter mass profile beyond i?out- The best fit parameters 
that minimize the of the comparison between the temper- 
atures predicted by Eq.Q and the observed temperature pro- 
file are searched in the intervals 50kpc < rg < 1.5 Mpc and 
2 < c < 13 for both SCDM50 and ACDM70 cosmologies. The 
results of the procedure are listed in Tabled We have also per- 
formed the same fitting procedure by neglecting the gas mass 
contribution to the total mass and find that in this case the fitting 
procedure yields substantially larger x^. We do not obtain any 
solution for the Virgo cluster, MKW 9, Hydra A and A 399, 
which are therefore excluded from the sample. For Virgo we 
obtain an unphysically small scale radius (and large concen- 
tration), because of the small field of view. The failure of the 
fit for the other clusters can be attributed to the low degree of 
smoothness of their temperature profiles (no smoothing of the 



data is used in our procedure). 

From the best fit parameters we compute ta using 

^ ^ 3 A/tot (< ta) 
47rpc,z?'A 



(3) 



where pc,z — {^H^) / {^nC) is the critical den- 
sity. We estimate ta for various overdensities: 
A = 2500,2000,1500,1000,500,200 for both cosmolo- 
gies, and at A = Ac, defined such that ta^ corresponds 
to the cluster virial radius rvii. For a SCDM cosmology, 
Ac = 178, while Ac = 178 + 82a; - 39a;2 for ACDM, where 
X = n{z) - 1 and n{z) = n^{i + z)V(f^m(l + z)^ + flA) 
( [Bryan & Norman (1998)| . In the computation of ta we use 
either 1) Mtot = ATdm + A/gas or, neglecting the gas mass 
contribution, 2) Aftot — A/dm- Assumption 1) naturally leads 
to values for r/\ that are larger than those computed assuming 
2). For both cosmologies, the relative difference between 
these values is small: 4-8 percent for A = 2500 and 10-15 
percent at A = 500, depending on the gas mass fraction of 
the cluster Even though assumption 1) is consistent with our 
procedure for the dark matter profile estimation, we use the 
TA values computed considering dark matter only for the 
following reason. When gas is taken into account, one needs 
to extrapolate Afgas(< r) for r > i?out when needed. We 
extrapolate using a least squares fit with an integrated /3-profile 
with the parameter f3 fixed to values taken from the literature: 



for Perseus (A 426) from Ettori et al. (1998) ] for A 311 2 from 
[Sanderson et al. (2003)[ and from .Peterson et al. (2003)[ for the 
remaining clusters. The mean relative error measured in the 
observed region is propagated to the extrapolated values. We 
find that for most of the clusters in our sample extrapolation 
is needed even for A — 1500 and that therefore the procedure 
is very unreliable at large radii. In addition, the extrapolated 
values for the gas mass are very sensitive to the value of 
the parameter f3. We therefore prefer, as widely done in the 
analysis of numerical simulations and in studies similar to our, 
assumption 2), i.e. we neglect the gas mass contribution when 
computing r a . 

The size-temperature relation ta oc < > predicted by 
self-similarity (Mohr & Evrard (1997)^ and its normalization 
provided by catalogues of simulated clusters allows an estimate 
of TA from the mean emission-weighted temperature < Tx > 
alone. We employ the widely used predictions given by the 
simulation of [Evrard et al. (1996)[ and compute the characteris- 
tic radii r^. We estimate the mean emission-weighted cluster 
temperature < Tx > by excluding the cooling flow region (i.e. 
we excise bins with cooling time less or equal to 15 Gyr) and 
fitting the remaining emission-weighted temperature bins with 
a constant. We list the mean emission-weighted temperatures 
< Tx > in Tableia For both SCDM50 and ACDM70, we find 
that is systematically larger than the values we compute 
by directly using the best fit NFW profiles. This is independent 
on whether our estimate of ta is smaller or larger than i?out 
and we find that the relative difference can be as high as 60-80 
percent for A — 200. This result has also been found by 
[Sanderson et al. (2003)[ with a different method to determine 
the total gravitational mass. 

From the set of (c, rg) parameters acceptable at 1 a we com- 
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pute, for each radius ta, the maximum and minimum value 
for the total mass within ta- In Tables|2]and|3]we list ta, dark 
matter mass and gas mass within ta for A — 2500 and 500 for 
SCDM50 and ACDM70, respectively. The quoted error related 
to the mass estimate is the root mean square of the difference 
between maximum and mean, and the difference between 
mean and minimum. When extrapolation is needed, the values 
for the gas mass are computed using the least squares fit 
described above. We compared our estimates with results 
from the literature. For clusters with published quantities 
comparable to ours, we generally find good agreement: 

- A 262: [Finoguenov et al. (2001a)| using temperature pro- 
files derived with ASCA, estimate r^^Q = 860 ±170 kpc 
and Mt,ot(< rgoo) = 9.3 ± 2.7 x lO^^^^^ SCDM50). 
Our estimate of r^QQ (see Table |2j is lower probably due 
to the fact that contrary topinoguenov et al. (2001a)| we do 
not include the gas mass contribution in its computation. 
Consequently also Aftot(< fsoo) is lower, but when com- 
puting Mtot(< 860 kpc) = 6.7 ± 1.8 x lO^M© we find 
consistency for the total mass estimate. 

- 2A 0335H-096: as for A 262 the estimate for r5oo of 
p^inoguen ov et al. (2001a)| is at variance with our, but agree- 
ment is found when the total mass is compared at the same 
radiu s: 34.4 ± 11.3 x IO^Mq and 26.1 ± 3.64 x lO^^M© 
being lFinoguenove t al. (20 01a)[s and our estimate, respec- 
tively, at 1340 kpc. Ett ori et al. (2002b)| derived total and 
gas masses from BeppoSAX data at A = 2500 and 1000 
using a method similar to our, with the only difference be- 
ing that their NFW profile estimates dark matter and gas 
masses altogether and not dark matter only as in our analy- 
sis. Due to this difference and to the larger radial range they 
access (i?out = 917 kpc with respect to i?out = 510 kpc 
in our work), the best fit scale radius is larger (rg — 626 
kpc) than our estimate and consequently the concentration 
parameter smaller. The characteristic radii ta estimated by 
[Ettori et al. (2002b) are also larger than ours, again due to 
the fact that we neglect the gas mass in the computation. 
Nevertheless, we find good agreement between the total 
and gas masses when evaluated at the same physical radius 
rather than at the same overdensity. 

- MKW 3s: we find excellent agreement between the esti- 
mates of |Finoguenov et al. (200 la)^ rsoo — 1210±140 kpc 
and Mtot(< ^500) = 25.7 ± 4.3 x IO^^Mq, and ours (see 
Table|2]for our values of Mdm and Afgas). 

- A 2052: we find excellent agreement between the estimates 
of [Finoguenov et al. (200l"a)| r^oo = 1280 ± 280 kpc and 
Mtot (< rsoo) = 30.2 ± 9.6 X lO^^A/© and ours (see Table 

- A 4059: the values given by [Finoguenov et al. (2001a)| 
(rsoo = 1390 ± 300 kpc and Aftot(< rsoo) = 38.7 ± 
12.0 X IO^^Mq) are in agreement with ours (see TableO. 

- A 496: as for A 262 and 2A 03 35-H096 the estimate of 
rsoo of [Finoguenov et al. (20 01ay is at variance with our. 
Disagreement is still found for the total mass estimated for 
rsoo = 1480 kpc: 47.0 ± 6.7 x 10"A^q and 19.9 ± 9.7 x 
10^^ Mq are [Finogueno v et al. (20 01a7s and our estimate, 
respectively. As for 2A 0335H-096 our estimates can also be 



compared to [Ettori et al. (2002b)[ s results. We find discrep- 
ancy as in the case of 2A 0335H-096: due to the much larger 
radial range they access (i?out = 845 kpc with respect to 
^out = 319 kpc in our work), the best fit NFW scale radius 
is much larger (rg — 738 kpc) than our estimate (see Table 
Q and consequently the concentration parameter smaller. 
Furthermore, even when compared at the same radius, our 
results for both total and gas masses are w 50% smaller 
than [Ettori et al. (2002b)[ s. The disagreement might be due 
to the fact that the comparison is made at radii much larger 
than i?out, where the mass extrapolation is large. 

- A 3 1 12: we find excellent agreement between the estimates 
of [Finoguenov et al. (2001a)| rsoo — 1350 ± 120 kpc and 
Mtot{< T-soo) = 35.4 ±4.5 x lO^^Afg and ours (see Table 

la. 

- A 1795: agreement is found between the estimates of 
[Finoguenov et al. (2001a)[ rsoo = 2000 ± 560 kpc and 
Mtot(< rsoo) = 116.6 ± 47.2 x 10"A/q and ours (see 
Table|2j. Better agreement is found when masses are eval- 
uated at the same radius: we estimate Aftot(< 2000 kpc) = 
82.8 ± 10.2 X lO^^Afg. As for 2A 0335H-096 and A 496 
we find a difference between our best fit NFW parame- 
ters and those computed by [Ettori et al. (2002b)| (they can 
use i?out = 1584 kpc with respect to i?out = 900 kpc in 
this work). The difference between [Ettori et al. (2002b)| s 
characteristic radii and total and gas masses is neverthe- 
less small and good agreement is found when these quan- 
tities are computed at the same radii, i.e. Mtot{< 771 
kpc) = 40.0 ± 1.9 X lO^Af p and Afgas(< 771 k pc) = 
53.7 ± 0.7 X from Ett ori et al. (2002b)[ versus 
Aftot(< 771 kpc) = 37.5±2.5xlOi3Af0 andAfgas(< 771 
kpc) = 47.7 ± 0.8 X 10^2^/0 in this work. 

- A 426: even though there is a difference between the NFW 
parameters estimated by ,Ettori et al. (2002b), and ours we 
find excellent agreement between the values of the charac- 
teristic radii and both total and gas masses. 

- A 1835: .Allen etal. (200 1)[ derived total masses from 
Chandra data at A = 2500 for both SCDM50 and 
ACDM70 and using a method similar to our. For both 
cosmologies, we find good agreement with their best fit 
NFW profile parameters and estimated r25oo- The small 
differences, in particular the smaller values we obtain for 
^2500, are very likely due to the fact that in their method 
the gas mass contribution is also included in the compu- 
tation of r25oo- When total masses are compared at the 
same radii, good agreement is found, i.e., for SCDM50, 
Aftot(< 720 kpc) = 5.41+^:^^ xlO^^ A^o and A^tot(< 72 
kpc) = 4.56 ± 0.35 x IO^'^Mq being [Allen et al. (2001)[ s 
value and our estimate, respectively. 

- ForNGC 533, A 1837 and Sersic 159—3 a fair comparison 
with previous work is not possible. 

To reiterate, the remainder of our analysis of our final sam- 
ple of 13 objects relies on the validity of extrapolating pro- 
files fitted to the inner region of the data: for both SCDM50 
and ACDM70, the number of objects for which the outer ra- 
dius i?out is larger than (or equal to) r a at the 95 percent level 
of confidence is: 11 (NGC 533, A 262, A 1837, Sersic 159-3, 
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2A 0335+096, MKW 3s, A 2052, A 4059, A 3112, A 1795 and 
A 1835) for A = 2500, 7 (NGC 533, A 1837, Sersic 159-3, 
MKW 3s, A 4059, A 1795 and A 1835) for A = 1500, only 3 
(A 1837, Sersic 159-3 and A 1835) for A 1000 and none for 
A = 500. 

4. Scaled temperature profiles 

In the following we investigate the structure of the tempera- 
ture profiles for our final sample of 13 cooling flow clusters. 
In Figs.[T](for SCDM50) and|2|(for ACDM70), we present the 
deprojected radial profiles plotted against the radius in units of 
Tvir (~ riso for SCDM and w rioi for ACDM at z = 0), where 
the temperature has been normalized by the mean emission- 
weighted temperature < Tx >■ From a visual inspection it is 
evident that a temperature gradient is present at large radii and 
that when normalized and scaled by the virial radius, tempera- 
ture profiles are remarkably similar In addition almost all the 
individual profiles clearly show a break radius 7-bi , a decrease 
of temperature from rbr inwards and a decline at radii larger 
than rbr- In tiie following, we investigate the shape in detail. 

4. 1. The break radius 

In the following we compute the break radius of the individual 
clusters in units of rvir. For each cluster we divide the scaled 
temperature profiles in two radial intervals: from the innermost 
bin to the bin i and from the bin i + 1 to the outermost one. 
Temperature profiles in each of the two intervals are then fitted 
using straight lines, power laws and exponential functions. All 
the nine combinations are used and for each pair of fitting func- 
tions, i is varied until the best fit is achieved. Clusters that do 
not show a clear temperature decrease in the outer region are 
excluded. These are A 262, A 496 and Perseus. For the remain- 
ing clusters the bin i which gives the best fit is independent 
of the choice of the fitting functions. The break radius a;bi — 
Thi/fvir is then defined as {xi + Xi+i)/2 (where Xi is the dis- 
tance of the bin i from the center in units of rvir) and its uncer- 
tainty (xi+i — Xi)/2. InFig.|3lwe show the break radius x^°'' — 
^br°V^vir of the projected temperature profiles in SCDM50 
as a function of redshift. The mean value of x^°'^ is 0.11 
with a standard deviation of 0.01. Using a different method, 
jPe Grandi & Molendi (2002)] find 0.20 for flie mean value of 
^br°^ for a sample of 11 cooling flow and 10 non-cooling flow 
clusters and a lower value, 0.16, for the cooling flow clusters 
only. Taking into account fliat jPe Grandi & Molendi (2002)| 
who also used projected profiles, computed the break radius 
for cooling flow clusters by excluding the cooling region and 
fitting the profiles with a constant temperature for r < r^^"' 
and with a line for r > r?^°K the fact that their estimate is 

br ' 

larger than ours is not surprising. In ACDM70 the mean value 
of x^'^°' is reduced to 0.08 with a standard deviation of 0.01. 
In agreement with De Gra ndi & Molendi (2002)| we find that 
the intrinsic dispersion of the parent population of scaled break 
radii (assumed to be distributed like a Gaussian) is consistent 
with (for both SCDM50 and ACDM70). By performing the 
same analysis using deprojected temperature profiles we find 



0.12 and 0.09 for the mean value of x'^''^'°^ for SCDM50 and 
ACDM70, respectively, and standard deviation of 0.01. This 
latter results shows consistency between break radii of pro- 
jected and deprojected profiles. 

4.2. The outer region 

For comparison with other studies, we quantify the decline 
seen in Fig. ^ (i.e. in a SCDM50 cosmology) for radii larger 
than ~ 0.15rvit by simply modeling the data with a line. From 
the data shown in Fig. ^ we select, for each object, bins at 
distances larger than rji^°^™^ — 0.12 r^ir. We fit the selected 
points to a straight line and find a slope equal to— 0.97±0.14 
(X^ = 20.38, d.o.f. = 34). 

Since many authors use projected temperature profiles, we 
apply the same procedure to the projected temperature pro- 
files. These are scaled by riso and normalized by the mean 
emission-weighted temperature < r^'°-' >, i.e. computed 
using the projected data. In this case we obtain —0.94 ±0.11 

(X^ = 26.24, d.o.f. = 35) f or the slope. 

For a fair comparison with jPe Grandi & Molendi (2002)| we 
compute the value for the slope using their procedure: we use 
projected temperature profiles (normalized by < T^'^°' > and 
scaled by rf^g), exclude temperature bins within the cooling 
region and fit the data with the broken line model given in 
Eq. (1) in jPe Grandi & Molendi (2002)| i.e. with a constant ti, 
up to a break radius Xb and with a Une with slope m beyond 
Xf,. In excellent agreement with jPe Grandi & Molendi (2002)| 
we find h = 1.02 ± 0.02, x^ = 0.10 ± 0.05 and a slope 
m = —1.15 ± 0.10. We also model the profiles scaled in 
the SCPM50 cosmology with a polytropic equation of state, 
T/ < Tx To(l + (r/ax)2)-3/3(7-i)/2, take values 
for (3 and from the literature: for Perseus (A 426) from 
lEttorietal. (1W5]1 for A 3112 from [Sanderson et al. (2003)| 
and froni^ Petereon et al. (2003)| for the remaining clusters. We 
then use the mean values /3 = 0.69, Ox = 0.087 rigg in the fit. 
Forr > r^°/'°' we find Tq = 1. 36±0.07and 7 = 1.24± 0.03, 
in agreement with the results of pVIarkevitch et al. (1998)| 
Recent numerical simulation are in general able to produce 
temperature profiles quite similar to the observed ones at 
large radii (see jPe Grandi & Molendi (2002)' in particular 
Fig. 12, for a comparison with simulations results published 
before 2001). Recently, jPe Grandi & Molendi (2002)) s 
profiles have been compared with simulated pro- 
files of | Valdarnini {IbW^ [borgani et al. and 
jBorgani et al. (2004)| and for radii larger that « 0.15rvir 
it was found good agreement. Given the agreement between 
Pe Grandi & Molendi (2002) s profiles and ours at large 
cluster-centric distances, these simulations should provide 
a good description of the temperature profiles presented 
here as well. Here, we compare the projected temperature 
profiles with the numerical pr edictions by | Loken et al. (2002)| 
and find good agreement. jLoken et''al^2002)| compared 
their results and found good agreement with observations 
of jMarkevitchetal. (1998)1 and flie Be ppoSAX profiles pre- 
sentedliTn^nran^^^^ole ndi (20(3^ The comparison with 
|Loken et al. (2002), s results is shown in Fig.|4] where we plot 
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Table 1. The radius at the outer end of the last radial bin considered, i?out and the best fit parameters for the NFW dark matter 
density profile for SCDMH50 and ACDMH70. Virgo, Hydra A and A 399 do not converge in the rg — c plane. The quoted errors 
on and c are RMS of the 68 percent joint confidence limits. 
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Table!. Mean emission-weighted temperature < Tx > and characteristic radii ta, dark matter mass A/dm(< ''a) and gas mass 
Mgas(< r-A) for overdensities A = 2500 and 500 in the SCDM50 cosmology. 
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< Tx > 
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(keV) 
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NGC 533 


1.21(0.05) 
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0.81(0.04) 


0.56(0.08) 
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1.30(0.09) 


1.89(0.27) 


A 262 


2.17(0.01) 


322(21) 
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663(51) 


4.44(0.80) 


9.24(1.46) 


A 1837 


4.00(0.13) 
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9.02(1.77) 


9.03(2.16) 


1021(208) 


18.99(8.37) 


29.08(9.50) 


Sersic 159-3 


2.29(0.07) 


413(9) 


6.05(0.25) 


10.53(0.60) 


888(28) 


12.02(0.86) 


25.17(1.97) 


2A 0335+096 


2.94(0.05) 


439(12) 


6.79(0.35) 


11.69(1.13) 


947(32) 


13.63(1.03) 


26.02(2.79) 


MKW 3s 


3.48(0.14) 


503(11) 


10.59(0.44) 


13.59(0.76) 


1095(34) 


21.80(1.51) 


35.28(2.27) 


A 2052 


2.99(0.08) 


481(21) 


9.00(0.74) 


11.33(1.00) 


1068(65) 


19.63(2.60) 


31.64(3.15) 


A 4059 


3.96(0.05) 


640(75) 


21.88(3.77) 


18.30(3.41) 


1637(372) 


73.05(30.38) 


49.77(12.66) 


A 496 


4.20(0.19) 


420(48) 


5.93(1.37) 


9.65(1.96) 


880(123) 


10.88(3.55) 


25.39(5.53) 


A3112 


4.33(0.11) 


578(104) 


17.47(5.18) 


20.98(5.15) 


1321(390) 


41.64(25.59) 


54.49(20.59) 


A 1795 


5.62(0.04) 


631(16) 


22.02(1.02) 


36.29(2.09) 


1443(54) 


52.58(4.15) 


96.49(6.31) 
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6.78(0.37) 
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70.93(7.26) 
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the mean error- weighted projected profile for r > O.lrvir for 
our sample of 13 cooling flow clusters on top of the results 
shown in Fig. 6 of Loke n et al. (2002)| |Loken et al. (2002)1 
find that the simulated clusters in their catalogs (both ACDM 
and SCDM catalogs) show temperature profiles which can 
be well described by a universal temperature profile of the 
form T/ < Tx To(l + r/a^)-\ with To = 1.33, 

Ox = f'vii/l-S and S = 1.6. This profile is clearly not suitable 
to reproduce the profile typical of cooling flow cluster over 
the entire observed radial range and a central region must 
be excised. The radial range used by |Loken et al. (2002)| 
in the fit is (0.04 - 1.00) rvir. Due to the limited radial 
range of our data, we fit the universal profile found by 
[Loken et al. (2002)| to the projected temperature profiles for 
< rP™Vri8o >= 0.11 < r/nsQ < w/riso = 0.49 by 
assuming ax = ri8o/1.5. We find Tq = 1.31 ± 0.05 and 



5 = 1.67 ± 0.20. Despite the good agreement, we stress 
that a more robust comparison will be only possible with 
observations spanning a much larger radial range than the one 
used in the present work. 

4.3. The inner region 

The temperature decline towards the center is characteristic of 
cooling flow clusters. For the sample in this paper, the struc- 
ture of deprojected temperature profiles has been extensively 
analyzed in |Kaastra et al. (2004)| The parametrization: 

T(r)=T, + (T,-T,)^||^, (4) 

is found to provide a remarkably good description of the tem- 
perature decline, with Tc set equal to the temperature of the 



R. Piffaretti et al.: Temperature and entropy profiles of nearby cooling flow clusters observed with XMM-Newton 7 

Table 3. Characteristic radii ta, dark matter mass Mdm(< ?'a) and gas mass Mgas(< ta) for overdensities A = 2500 and 500 
in the ACDM70 cosmology. 
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Fig. 1. Scaled temperature profiles (deprojected) in SCDM50 cosmology: the radius is scaled with the virial radius r^\^ = rigo, 
while the temperature is normalized by the mean emission- weighted temperature < Tx >. Clusters are related to symbols as 
follows: NGC 533 {crosses), A 262 (Med squares), A 1837 (fii/ed diamonds), Sersic 159-3 {GUed circles), 2A 0335+096 (open 
triangles), MKW 3s (open pentagons), A 2052 (filled triangles), A 4059 (open diamonds), A 496 (open hexagon), A 3 112 (open 
stars), A 1795 (open squares), Perseus (open circles) and A 1835 (fflJed stars). 



irmermost bin and /i = 2. In SCDM50 a correlation be- 
tween the fitting parameters Th and Tc of the form Tc = 
(4.1 ± 0.8)r^^-^''^° ", with rc in kpc and Th in keV, is found. 

1/2 

Assuming that ta scales with ' one concludes that the char- 



acteristic radius Tc does not simply scale with the virial radius. 
In the following we further explore this implication by consid- 
ering the mean emission-weighted temperature < Tx > in- 
stead is Th and using the estimated ta instead of assuming 
r-A oc tI''^. For both SCDM50 and ACDM70 cosmologies we 
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find that the characteristic length scale Vc is not linked to the 
size of the dark matted halo hosting the cooling flow in a simple 
way. In SCDM50, for example, we find rc (x< Tx >i-78±o.io 
and r25oo oc< Tx >"-88±o.o2^ ^^^^ ^^^^^ j,p^ 

< Tx > in keV. In this context we notice that the flatten- 
ing of the composite scaled temperature profile in the range 
~ 0.1 - 0.2rvir in SCDM50 (see Fig.[T) is artificial, as it is due 
to the fact that temperature profiles in the inner region do not 
simply scale with the virial radius. 

IKaastra et al. (2004)| searched for correlations between param- 
eters of the cooUng gas without finding any significant one 
l |Kaastra et al. (2004)| section 4.1). Having additional informa- 
tion about the characteristic radii ta, we investigated correla- 
tions of these with the parameters of the cooling gas, such as 
relative temperature decrement and cooling radius. We do not 
find any significant correlation. 

As shown in the previous section, for radii larger than ^ 
0.15rvir most of the numerical simulations are able to produce 
temperature profiles quite similar to the ones presented in this 
work. The fair comparison with simulated cluster temperature 
profiles in the central region, i.e. for radii smaller than the break 
radius rbr, is not possible yet. In fact, the central temperature 
drop has not been reproduced consistently by any simulation 
(e.g., ITornatore et al. ( 2003) Borgan i et al. (20 04)' and refer- 
ences therein). In this context, our temperature parametrization 
for the inner region of cooling flow clusters (Eq.|3 and the cor- 



relation between its free parameters (rc (x< Tx >i 'i'8±o iO) 
provides an important test for forthcoming simulation models. 



5. Entropy profiles and the L-T relation 

The entropy is a useful quantity to probe the thermodynamic 
history of the hot gas in groups and clusters. As it has become 
customary, we define the entropy as 5 = kT/ric^^, where T 
and rio are the deprojected electron temperature and density, 
respectively. S, which is by some authors referred as to 
'entropy index', is related to the thermodynamic entropy via a 
logarithm and an additive constant and therefore conserved in 
any adiabatic process. 

Particular interest has been given to the shape of entropy 
profiles in groups, since preheating models predict both large 
isentropic cores and a break in the entropy temperature relation 
(e.g., |Tozzi & Norman (200 1)^ . While there is evidence that 
even at group scale no large isentropic cores are present 
( jPratt & Arnaud (2003)| jPonman et al. (2003)> , the entropy 
profiles of cooling flow clusters are expected to increase 
monotonically moving outwards and to show no isentropic 
cores (e.g., |McCarthy et al. (2004)| >. 

In Fig.|5]we plot entropy profiles for the 13 cooling flow clus- 
ters with the radial coordinate scaled by rvir in the SCDM50 
cosmology. A visual inspection of the individual profiles in 
the plot shows that, as expected for systems with short cooling 
times in their central regions, no cores are present and the 
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Fig. 3. Break radii rj^'.°^ in units of r^h- as a function of redshift in SCDM50. The horizontal line is the eiTor- weighted mean 
value. 



entropy of the ICM increases monotonically moving outwards. 
In the standard self-similar scenario, one predicts 
S ( X h-^/^{z)T, where h^{z) = + z)^ + Ha 

(e.g., Pratt & Arnaud (2004)). On the other hand, recent results 
by .Ponman et al. (2003) suggest that entropy scales with 
temperature as 5 oc T^-^^^ the so-called "entropy ramp", 
instead of S (x T. Both scalings S oc h^^/^{z)T and 
5* cx r"-^^ considerably reduce the dispersion of 

the profiles shown in Fig. |5l Using Tvir to scale radii, we 
quantify the dispersion of the scaled profiles using the standard 
deviation and the mean of the scaled entropy values at a 
fixed scaled radius. As mean cluster temperature T, the mean 
emission-weighted temperature < Tx > (see Table|3 is used. 
The scaled entropy values are evaluated at fractions of rvir for 
which no extrapolation is needed and we compute the ratio 
between their standard deviation a and mean m to quantify the 
relative dispersion of the scaled profiles. For both SCDM50 
and ACDM70 cosmologies we find that the dispersion is less 
if the "entropy ramp" scaling is used. For SCDM50 we find 
cr/m = 0.37 and ct/to = 0.29 at O.lrvh- for the S" cx T and 
S cx T"^^ temperature scaling, respectively. The relative 
difference of the cr/m ratios is the same in ACDM70 and 
the relative dispersion cr/m varies with radius only by a few 
percent (e.g. 5% between 0.03 and 0.2rvir in SCDM50). 
From scalings of the entropy profiles of 5 clusters (A 1983, 
A 2717, MKW9, A 1991 and A 1413), |Pratt & Arnaud (2004)| 
also find that the scatter in the scaled entropy profiles is 



smaller if one uses'Pon man et al. ( 2003)' s scaling rather than 
the standard self-similar one. In addition to the dispersion 
analysis, we attempted to gain further information on the 
scaling of entropy with temperature by directly fitting the 
S-T data (i.e. 5(0.1 r2oo) versus < Tx >), but due to the 
small sample size and the fact that the systems do not span 
a wide temperature range, we do not obtain any statistically 
significant constrain on its slope. The profiles scaled using 
the relation S cx h-'^'^{z)T°-^^ are shown in Fig. |6l for 
SCDM50. We fit the data shown in Fig.|6lwith a line in log-log 
space (with errors in both coordinates) and find, in excellent 
agreement with the value found by |Pratt & Arnaud (2004)| 
a slope equal to 0.95 ± 0.02. Similarly, using Chandra data, 
[Ettori et al. (2002a)| find a slope equal to 0.97 for the entropy 
profile of A 1795. Therefore, our analysis gives additional evi- 
dence for a slope close to, but slightly shallower than 1.1, the 
value predicted by shock dominated spherical collapse models 
( |Tozzi & Norman (2001)| l. This behavior is also a common 
feature of numerical simulations (see |Pratt & Arnaud (2004)| 
and references therein). 

From the scaled entropy profiles we can obtain an estimate 
of the normalization of the S cx T°^^ relation found by 
[Ponman et al. (2003)| For the entropy profiles with entropy 
values scaled as 5 cx h~'^/^(z) [T /IGVeVf-'^^ and the 
radial coordinate scaled by r2oo we find a mean value at 
0.1 r2oo equal to 504 (with standard deviation 140). The latter 
value is in excellent agreement with the normalization of 




0.2 0.4 0.6 0.8 1.0 1.2 

Fig. 4. Mean error-weighted projected temperature profile for our sample (filled triangles, with 1 a errors) as a function of radius 
in units of r^ir together with the results shown in Fig. 6 of Loken et al. (2002). The two solid lines define the approximate band 
that encloses most temperature profiles and their error bars of Markevitch et al. (1998). 



IPonman et al. (2003)| s 5(0.1 raoo)-^ relation at T = lOkeV 
(see their Figure 5). 

Although the relation S oc T'^-^^ provides evidence for 
similarity-breaking, the most robust evidence for the lack 
of self-similarity in the ICM properties is provided by 
the luminosity-temperature (L-T) relation. The L oc T" 
law, with a = 2 expected for self-similar systems is 
at variance with observations at both cluster and group 



scales. For clusters, most authors (e.g. [White et al. (1997)| 
[Arnau d & Evrard (1999)1 have found a ^ 3, while for groups 
there is evidence for a much steeper slope a ^ 5 (e.g. 
elsdon & Ponman (2000a}| |Relsdon & Ponman (2000E] 
ue & Wu (2000)1. The latter result is however not confirmed 
by the work of ,Mulchaey & Zabl udoff (1998)' and most im- 
portant by the recent analysis by 'Osmond & Ponman (2004)| 
as both works indicate that the L-T relation at group scales 
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Fig. 5. Entropy profiles in SCDM50: the radius is scaled with the radius Tvir = ''iso^ while no scaling has been applied to the 
entropy values. Clusters are related to symbols as follows: NGC 533 (crosses), A 262 {filled squares), A 1837 (filled diamonds), 
Sersic 159-3 (filled circles), 2A 0335+096 (open triangles), MKW 3s (open pentagons), A 2052 (filled triangles), A 4059 (open 
diamonds), A 496 (open hexagon), A 31 12 (open stars), A 1795 (open squares), Perseus (open circles) and A 1835 (filled stars). 
The dashed line indicates the best fit power-law with the power index value fixed to 1.1 predicted by |Tozzi & Norman (200 1)| 



is consistent with and extrapolation of the trend observed at 
cluster scale ( a ~ 3). 

A great deal of semi- analytical and numerical studies 
have attempted to give an explanation to the observed 
similarity-breaking. While it is clear that models incorpo- 
rating gravitationally-driven processes only fail to match 
observations, there is evidence that models incorporating a 
mixture of both cooling and heating processes are likely to be 
successful (see Ponm an et al. (20 03) and references therein). 
Such models should of course also explain the properties of 
cooling flow clusters. Cooling flow clusters are known to 
enhance the scattering of the observed L-T relation, which 
is reduced once the contribution from the cooling region is 
excised or cooling flow clusters are excluded from the sample 
(Markevitch (1998) Ettori et al. (2002b)i. On the other hand, 
the scatter in not random and theoretical models incorporating 
gas cooling should be compared to "uncorrected" data. This 
kind of comparison has never been presented in semi-analytical 
or numerical works until recently in [McCarthy et al. (2004)| s 
work. [McCa rthy et al. (2004)| have extended the work by 
[Babul et al. (2002)] to include radiative cooling and compared 
their predictions with observations without any "cooling flow" 
correction. 

In Fig. |7] we show the L-T relation in SCDM50 for the 



13 cooling flow clusters we have analyzed. According to 
[McCarthy et al. (2004)[ s predictions, clusters with large 
mass deposition rates M lie on the high luminosity side 
of the scaling relation. We take mass deposition rates from 
[Peterson et al. (2003)[ and, for the objects not present in the 
latter work (A 3112 and A 426) from iPeres et a l.(1998)[ 
Mass deposition rates for A 1837 are not available in the 
literature. We divide the clusters into three groups: those with 
M > SOOMQyr^^ ("massive cooling flow clusters"), with 
M > lOOMQyr^^ ("moderate cooling flow clusters") and 
clusters with mass deposition rates consistent with zero. Even 
though our sample is small, a clear trend consistent with the 
findings of pVIcCarthy et al. (2 004)' can be seen from a visual 
inspection of Fig. [71 in which we also show the predictions 
by [McCarthy et al. (2004)[ for an entropy injection level of 
200 keV cm"^ and 500 keV cm^ (models evolved for 10 Gyrs, 
see [Babul et al. (2003) ) . "Massive cooling flow clusters" are 
consistent with being systems that experienced "mild" (~ 200 
keV cm^) levels of entropy injection. 

We feel it should be pointed out that the entropy injection lev- 
els in McCarth y et al. (2004)[ s model are quite high compared 
to most of the preheating models (e.g., see [Voit et al. (2003)[ 
and references therein). However, a comparison between 
preheating models is beyond the scope of our work, being our 
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Fig. 6. Same as Figure|5lbut with the scaling S oc h^'^/'^{z) T*'-^^ used instead. The solid line indicates the best fit power-law 
(best fit power index 0.95) and the dashed one the best fit power- law with the power index value fixed to 1.1. 



main purpose that of providing "uncorrected" L-T data for a 
sample of cooling flow clusters. In Fig.|8]we show a second 
relationship predicted by the work of McCarth y et al. (2004)| 
the correlation between core radius and central entropy. As 
expected, we find that the entropy at the core radius ax 
increases with increasing core radius, with the exception of the 
low mass system NGC 533. Using core radii taken from the 
literature (see above) and excluding NGC 533 from the fit, we 
find: log[S'(ax)] = 1.481og[ax] - 0.91 (entropy in keV cm^ 
and core radius in kpc). 

6. Summary and Conclusions 

We have investigated temperature and entropy profiles of cool- 
ing flow clusters using the results from the spatially resolved 
spectra taken with the EPIC cameras of XMM-Newton pre- 
sented in Ka astra et al. (2004)| 

Using deprojected plasma densities and temperatures we com- 
pute the total gravitational mass profile and the virial radius of 
each clusters and find that: 

- The inclusion of the mass contribution from the gas in the 
estimate of the total mass of the cluster is very important 
and gives more robust estimates. 

- As recently found by other authors, the estimated virial 
radii substantially deviate from those estimated assuming 
the standard size-temperature relation. 



- Our total and gas mass estimates are in good agreement 
with published ASCA and BeppoSax estimates. 

Our main findings concerning the temperature profiles are: 

- When normalized and scaled by the virial radius, the tem- 
perature profiles are, except in the central cooling region, 
remarkably similar. 

- The profiles show a break radius rbi ^ 0.12 and 0.09 Tvir 
in SCDM50 and ACDM70 cosmology, respectively. 

- Most important, the temperature profiles show a decUne be- 
yond the break radius. 

- Our findings for the shape and slope of the profiles in this 
range is in excellent agreement with previous studies by 
IMarkevitch et al. (1998). and ,De"&andi & Molendi (2002)| 
In SCDM50, a linear fit to the deprojected temperature pro- 
files with cluster-centric distance scaled by the virial ra- 
dius, gives a slope of ~0.97± 0.14 between ~ O.lSrvir and 
0.50r™.. 

- For radii larger than ~ 0.15rvii (in SCDM50) most of the 
numerical simulations are able to produce temperature pro- 
files quite similar to the observed ones. 

- The temperature decline from rbr inwards is very well 
parametrized using a profile with a core (see Eq.^J. We 
find that Tc oc< Tx >i 'i'8±o io^ vvith the mean cluster tem- 
perature < Tx > in keV and Vc in kpc, and that Vc is not 
proportional to the cluster virial radius. 

We have studied entropy profiles and find that: 
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Fig. 7. The cooling flow "uncorrected" L-T relation in SCDM50. The plotted quantities are the mean emission-weighted tem- 
perature < Tx > and the bolometric X-ray luminosity within rspo- Filled triangles and open boxes indicate clusters with 
M > 300MQyr^^ and M > lOOMQyr"^, respectively. Mass deposition rates for A 1837 are not available in the litera- 
ture, while M for A 262 and NGC 533 are consistent with zero. The dashed line shows jMcCarthy et al. (2004), predictions for an 
entropy injection level of 200 keV cm^, while the dot-dashed line an entropy injection level of 500 keV cm^. 



- As expected for cooling flow clusters, i.e. massive systems 
with strong cooling, no isentropic cores are present. 

- The entropy of the gas increases monotonically moving 
outwards almost proportional to the radius. 

- The scatter in the scaled entropy profiles is smaller if the 
empirical relation S oc T°^^ jPonman et al. (2003)> is 
used instead of the standard self-similar relation S ocT. 

- The profiles are close to, but slightly shallower (ex r°-^^ 
with entropies scaled using S oc T*^^^) than those predicted 
by analytical models of shock dominated spherical collapse 
(cx ri i). 

- The normalization of the S-T relation derived from the 
scaling of the entropy profiles is in excellent agreement 
with the one found by Ponman et al. (2003)| 

- We confirm the existence of a tight correlation between 
core radius and central entropy. 
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